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The spatial correlation with classical lights, which has some similar aspects as that with entangled 
lights, is an interesting and fundamentally important topic. But the features of high-order spatial 
correlation with classical lights are not well known, and the types of high-order correlations pro¬ 
duced are of limit. Here, we propose a scheme to produce third-order spatial correlated states by 
modulating the phases of three laser beams. With the scheme we can produce Greenberger-Horne- 
Zeilinger-type (GHZ-type) and W-type spatial correlations with different phase modulations. Our 
scheme can be easily generalized to produce Worder spatial correlation states and to probe the 
aspects of different multi-partite spatial correlations. 
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The high-order spatial correlation effect has been of 
great importance in optics [l| . The investigations on this 
effect not only can he^ us to understand the nature of 
light more effectively but also make it advantageous 
to bring such effect into practice [ikl. 

In 1995, second-order spatial correlation phenomenon 
was observed experimentally with entangled photon pairs 
from spontaneous parametric down-conversion (SPDC) 
Hi , which led to the researches about the theories and 
applications of the second-order spatial correlation effect 
with entanglement IMS- Several years later, it was dis¬ 
covered that many of the features of the second-order 
spatial correlation obtained with entangled photon pairs, 
such as ghost imaging, subwavelength interference and so 
on, can also be realized with classical lights likl. Then, 
the third-order spatial correlation effect has been exten¬ 
sively studied with both entangled and classical lights 
|15l-[^ . But the features and types of high-order spatial 
correlation are not well known to us for both entangled 
and classical lights. As we know, for entangled lights, 
there are GHZ and W entangled states in three-partite 
system and more complex types of entangled states in 
multi-partite system [^, . The different types of en¬ 

tangled states can produce corresponding spatial corre¬ 
lations mililli- For classical lights, we do not know 
how to produce GHZ-type and W-type correlations in 
the three-partite system, and do not know the types of 
correlation in general multi-partite systems. 

In this work, we focus on high-order spacial correla¬ 
tion with classical lights. First we define the GHZ-type 
correlated state and W-type correlated state. Then, we 
present a scheme to produce GHZ-type and W-type cor¬ 
related states with three classical light sources generated 
by applying three sets of time-variable phase masks onto 
three laser beams, respectively. It is shown that when the 
three sets of phase masks are different but correlated, the 
state of the three classical lights is a GHZ-type spatial 
correlated state. When the three sets of phase masks are 
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identical, the state of the three classical lights is a W-type 
spatial correlated state. This scheme can be generalized 
to multi-partite system, which means that we can pro¬ 
duce multi-partite correlated states, and achieve different 
higher-order spatial correlations. 

We define the GHZ-type correlated state as a state 
with which one can achieve third-order correlation, but 
no second-order correlations, and the W-type correlated 
state with which one can achieve both third-order and 
second-order correlations. It should be noted that the 
definition of GHZ-type (W-type) correlated states are 
just similar to the GHZ-type (W-type)entangled states in 
which there are three-partite entanglement but no (have) 
two-partite entanglement. Furthermore, for multi-partite 
systems the types of correlations are not known as the 
entanglement does. Although the features of correlated 
states may be simpler than those of the entangled states, 
they are still unclear to us and the correlated states have 
many theoretical and experimental applications, such as 
non-locality without entanglement and correlation imag¬ 
ing. 

With the definition above, we present a scheme to pro¬ 
duce high-order correlation. A schematic diagram of the 
model for describing the third-order spatial correlation 
effect is shown in Fig. [H Three classical light sources 
are involved. Each classical source is obtained by ap¬ 
plying a set of phase masks onto a laser beam. These 
phase masks can be achieved with a spatial light mod¬ 
ulator (SLM). The values in these phase masks k) 
for the n-th (n = 1,2,3) classical source are spatially 
random-distributed and time-varying, with ^ being the 
position vector on the phase masks, and k = 1, • • • 
being the index of the phase modulation samples. Here, 
we consider that the three light sources are placed at 
one same coordinate system. It will be shown that we 
can obtain GHZ-type and W-type spatial correlations 
among the three classical sources by properly controlling 
the three phase masks. 

Here, we consider the single photon model of each laser 
beam. Then, the state of the laser beams can be written 
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FIG. 1. (color online). Standard schematic of third-order 
spatial correlation using classical lights. Laser: laser light 
source; SLM: spatial light modulator; detector: the detector 
to measure the field pattern of the light source. 


as 

= I d^|l(”)( 0 ),(n = 1 , 2 , 3 ), ( 1 ) 

where \ presents the state that one photon is emitted 
from the position ^ in the n-th laser beam. The state 
in Eq. (HD means that the single-photon state in every 
position ^ is coherent. After the phase modulations, the 
three-photon density matrix to describe the three light 
sources can be written as [ 28 | 
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(exp{i [tpi{Ci,k) + ip2{C2,k) 


+ ^ 3 (^ 3 , k) - ip3{U,k) - ip2{C5,k) - ipi{C 6 ,k)]}), 


( 2 ) 


where (■) = ^ defined to be the sample av- 

erage over the K phase modulations. To investigate the 
second-order spatial correlation, two-photon state should 
be considered, which can be obtained by tracing away one 
photon of the three-photon state. Without loss of gener¬ 
ality, the third photon is traced out, and the two-photon 
density matrix describing the other two light sources is 


P2 



( 3 ) 


Furthermore, the fields at any distances dn from the three 
light sources can be obtained using the Fresnel diffraction 
integral, respectively mil: 

= J d^an{Ohn{Xn,C), ( 4 ) 

where an(C) is the annihilation operator in the n-th op¬ 
tical source. Moreover, hn{xn^C) is the Green’s function 
associated to the propagation of the field from the source 
to the detection plane 

hn{Xn, e oc exp [i7r(Xn - O'^/Xdn] , ( 5 ) 


where A is the wavelength of the laser light source, and 
Xn is the position vector at the detection plane. Then, 
by cross-correlating the field patterns, we can obtain the 
third-order and second-order spatial correlation functions 

Q: 

Gsixi ,X2,X3) =Tr[p3£;L^ {xi {x2)E^~^ (x^) 

X Ei+\x 3 )Ei+\x 2 )E[+\x^)], 
G2ixi,X2) =Tr[p2E[~\xi)E^~\x2) 

X Ei+\x2)E[+\x,)]. 

In this section, we will consider the situation of third- 
order GHZ-type spatial correlation. Here, the modula¬ 
tion method is considered: the values in the phase masks 
and are statistically independent, uni¬ 

formly and randomly distributed in the range [ 0 , 27 r), 
while keeping = ip2{^^k) + (ps{^,k). That is, 

the three sets of phase masks loaded on the three laser 
beams are different but correlated. In this case, we get 
obviously 

(exp[i(^„(^, k)]) = 0 ,{n = 1,2,3), 

{exp{i[(fin{^,k) - <Pn’{^',k)]}) = Sn,n’d{^ - C')> 
(exp{i[<pi(^i,fc) - (p2{C2,k) - ‘P 3 {C 3 ,k)]}) 


In this way, the density matrix of the three-photon state 
can be written as 

P3 OC J d^id^2d^3\i^3{^l,^2,^3)){i^3{^l,i2,^3)\‘ ( 8 ) 

Here, in the case that 7^ ^2 = fhe term 

I'03(Cf 6,C3)) can^ denoted by |' 03 (Cf 6 )), which is 
a GHZ-type state [ 25 |: 


|' 03 (Cf 6 )) — 


V2 




+ 





( 9 ) 


In the other cases, C2, C3)) = 


'g’ 


,( 2 ) 
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( 3 ) 


is a product state. The GHZ-type state in Eq. m re¬ 
sults in the superposition of two different but indistin¬ 
guishable alternatives: hi{xi,^i)h2{x2,^2)h3{x3,^2) and 
C2)^2(^2 ,^i)^ 3(^3 5Cl): f^LC so-callcd thrcc-photon 
interference, as shown in Fig. [ 2 j The third-order spa¬ 
tial correlation function corresponding to the GHZ-type 
state in Eq. is 


G3{xi,X 2,X3) (X j diidi2\hl{xi,ii)h2{x2,^2)h3{xz,i2) 

+ hi{xi,£,2)h2{x2,S,i)hz{xz,ii)\^. 


(10) 


When the distances di , d2 and ds satisfy the condition: 
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^ ^ = 0, 

d 2 as 


(11) 
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(a) 


(b) 


FIG. 2. The representation of the three-photon inter¬ 
ference. It is the result of the superposition between 

these two indistinguishable three-photon amplitudes: (a), 
^ 2 )fe(iC 3 , ^ 2 ): the probability amplitude that 
a photon in position of field 1 goes to detector Di, a pho¬ 
ton in position ^2 of field 2 and 3 goes to detector D 2 and D3, 
respectively; (b), the probabil¬ 

ity amplitude that a photon in position ^2 of field 1 goes to 
detector Di, a photon in position of field 2 and 3 goes to 
detector D 2 and D3, respectively. 


Eq. (pT|) can be simplified as 


G 3 (a;i,a; 2 ,a; 3 ) oc 1 + 
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( 12 ) 


where D is the size of the SLM. The term in Eq. dUD 
means a perfect spatial correlation among the light at the 
three detection planes (xi, X2 and X3), which has possible 
applications, such as ghost imaging. 

By tracing out the third photon in the three-photon 
density matrix in Eq. (jS]) , the two-photon density matrix 
describing two of the three light sources is obtained: 


P2 






( 13 ) 


This means that the two-photon density matrix is 
a product of the two single-photon density matrices 

f d^n ~ 1:2). As a result, there is no 


second-order spatial correlation. 

In order to verify the third-order GHZ-type spatial cor¬ 
relation effect, we performed numerical simulations, and 
the results are shown in Eig. [ 3 l As the input for sim¬ 
ulations, we considered the following case. Three laser 
light sources with wavelength A = 532 nm illuminate 
three SLMs with size D = 2 mm, respectively. Then, the 
beam from the first source transmits through a three- 
point object shown in Eig. [ 3 ]^a), and a single-pixel detec¬ 
tor collects the light transmitted from the object. The 
intensities of the other two light beams are measured by 
two scanning detectors. The distances are set as di = 10 
cm, d2 = ds = 20 cm, which guarantees the spatial cor¬ 
relation at the three detection planes. Erom Eig. [ 3 fb) 
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FIG. 3. (color online). Simulated results of third-order spatial 
correlation effect with an object of three points, (a) is the ob¬ 
ject to be imaged; (b) and (c) are the images reconstructed by 
measuring the third-order spatial correlation functions with 
the two types of measurement operations, respectively; (d) 
and (e) are the results obtained by measuring the second- 
order spatial correlation functions, respectively. 


and Fig. [ 3 ]^c), it is found that the information of the 
object can be obtained by the third-order spatial cor¬ 
relation functions. In addition, there are two kinds of 
measurement operations to reconstruct the image of the 
object. The first way is to move one detector (xs) with 
the other detector {x2 = 0) fixed, and the magnification 
factor of the image is ds/di = 2 with the image shown 
in Fig. [ 3 ]^b). The second way is to move both detectors 
{x2 = xs) together, with the result shown in Fig. [ 3 fc), 
and the magnification factor of the image is 1 . However, 
the information of the object can not be reconstructed by 
second-order correlation functions, as shown in Fig. [Sfd, 
e). The simulation results above show the GHZ-type spa¬ 
tial correlated state can be attained with the modulation 
method, which can achieve third-order GHZ-type spatial 
correlation. 

Now we will consider the W-type spatial correlation 
case where the three sets of phase masks loaded on 
the three SLMs satisfy the condition: the three sets of 
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phase masks are the same, namely k) = k) = 
^ 3 (^ 5 ^) — l^he values of the phases added 

on the pixels in the SLMs are statistically independent, 
uniformly distributed in the range [0,27r). This case is 
equivalent to the model of observing one type of third- 
order spatial correlation effect with classical light in Ref. 

As these phase masks obey the relationships: 
(exp[z(^(^,fc)]) = 0 and (exp{i[ 99 (^ 1 , fc) -<^ 9 ( 6 ,^)]}) = 
S{^i — ^ 2 ), the three-photon density matrix to describe 
the three light sources is simplified to 

P3 OC J d^ld^2d^3|'03(Cl:65C3))('03(Cl5 6 ^ 6 ) 1 - ( 14 ) 

Here, in the case that 7 ^ ^2 7 ^ ^ 3 ? the term 

|'03(^i5 C 25 C 3 )) can be taken as a W-type state [26|: 

te({..6,&))=7sn|4?,..). (15) 

p n=i 

where the sum runs over the all 3! possible permuta¬ 
tions P of the set of integers 1,2,3. In the case that two 
of the three parameters ^ 1 , ^2 and ^3 are the same and 
different from the other one, such as 7 ^ ^2 = ^ 3 ^ the 
term |'03(^i,^ 2 ,C 2 )) can also be taken as a W-type state 
0 : 


I^3(ei,e2,6))=- 



l(2)\ 

/ 

/ 


+ 



+ 




I (3) 

'g>> 


I (3) 


(16) 


In the other cases, the states |'03(Ci5 ^ 2 , C 3 )) = 

^[ 1 ^) ^[ 2 ^) ^[ 3 ^) product states. Therefore, the 


three-photon density matrix describing the three light 
sources is the incoherent superposition of a set of W-type 
and product states |' 03 (Ci 5 C 25 ^ 3 ))? leading to third-order 
spatial correlation. 

By tracing the third photon in the three-photon den¬ 
sity matrix in Eq. m, we obtain the two-photon density 
matrix: 


P2^ Idild^2\M^l,^2)){M^l,^2)\, (17) 

where the terms \ip 2 {^i,^ 2 )) in Eq- itTfll are Bell-type or 


product states 0 : 


|V’2(Cl56)) 


V2 


EpU 


2 

n=l 




1 


Cl 7^ C 2 , 

others. 


(18) 


This means the second-order spatial correlation still ex¬ 
ists in this case. Hence, the W-type spatial correlated 
state is produced, which can realize the W-type spatial 
correlation effect. 

In the analytical results above, it shows that the third- 
order GHZ-type and W-type spatial correlations can be 
realized with three classical light sources. Without loss of 
generality, this scheme can be generalized to achieve N- 
order GHZ-type and W-type spatial correlations with N 
classical light sources. If the N sets of phase masks satisfy 
the condition: k) = ^^=2 k) one can achieve 

the Worder GHZ-type spatial correlation. If the N sets 
of phase masks loaded on the N laser beams are the same 
one can achieve the Worder W-type spatial correlation. 
Specially interesting, if the N sets of phase masks sat¬ 
isfy other conditions, we may produce the Worder spa¬ 
tial correlation with the types more than GHZ-type and 
W-type spatial correlations. This means that we may 
produce the inequivalent multi-partite correlated states 
and different types of spatial correlations, and probe the 
features and the applications of multi-partite correlation 
with classical lights. 

In conclusion, we propose a scheme to produce the 
high-order spatial correlation with classical light sources. 
We show that GHZ-type and W-type correlated states 
can be produced by applying three sets of time-variable 
phase masks onto three laser beams, respectively. When 
the three sets of phase masks are different but correlated, 
the state of the three classical lights is a GHZ-type spa¬ 
tial correlated state. When the three sets of phase masks 
are identical, the state of the three classical lights is a 
W-type spatial correlated state. This scheme can be eas¬ 
ily generalized to multi-partite system to produce more 
types of multi-partite correlated states, and achieve more 
kinds of high-order spatial correlations besides the GHZ- 
type and W-type correlations. These results can deepen 
the basic understanding of the high-order spatial correla¬ 
tion effect effectively which is fundamentally important. 
Furthermore, it may lead to novel high-order correlation 
effects and applications. 
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